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Abstract 

The quantum dynamics of the symmetry broken A (<i>^)^ scalar field theory in the 
presence of an homogeneous external field is investigated in the large N limit. We 
choose as initial state the ground state for a constant external field J. The sign of 
the external field is suddenly flipped from J" to — J" at a given time and the subse- 
quent quantum dynamics calculated. Spinodal instabilities and parametric resonances 
produce large quantum fluctuations in the field components transverse to the external 
field. This allows the order parameter to turn around the maximum of the potential 
for intermediate times. Subsequently, the order parameter starts to oscillate near the 
global minimum for external field — J7, entering a novel quasi-periodic regime. 

PACS: ll.lO.-z, ll.15.Pg, ll.30.Qc 



Contents 



I Introduction 



2 


Classical field dynamics 






2.1 Classical evolution equation and initial conditions 




2.2 The classical dynamics 



1 



3 Quantum field dynamics 



3.1 Quantum equations of motion and initial conditions| 7 

3.2 Early time quantum evolution 9 

10 

10 



3.3 Intermediate time quantum evolution] 11 



3.4 Constraints on the trajectory due to the conservation of the energy] 14 

i Comments and Conclusions 



5 Acknowledgements] 



6 Appendix] 



15 



16 



16 



^.1 


Appendix A: Spinodal instability 




D.2 


Appendix B: Ground state results 





1 Introduction 

The understanding of the physics at large energy density situations is fundamental to achieve 
a proper picture both for ultrarelativistic heavy ion collisions [Q] and for the early universe |Q 
^, H . These physical phenomena cannot be treated with the usual quantum field perturbation 
methods and call for self-consistent nonperturbative methods as the large approach and 
Hartree approximations [0, ]^, |, |, ]^, |]. 

Here, we study effects of external fields in non-perturbative quantum field dynamics using 
the large N limit method. 

The effective potential of the 0{N) A theory with broken symmetry has two degener- 
ate minima for zero external field. In the presence of a small uniform external field the 
potential becomes tilted and the degeneracy between the minima lifted. The absolute min- 
imum becomes the vacuum while the other local minimum becomes a metastable vacuum. 
See Fig. |. 

We take as initial state the ground state. That is, the field expectation value is initially 
at the absolute minimum and the field modes are in the vacuum. Then, at a given time 
we change the sign of the external field: — J ■ Flipping the sign of the external field 

exchanges the local minima. That is, the absolute minimum becomes the false vacuum and 
viceversa. As a consequence, the system is now near the false vacuum. Spinodal resonances 
develop and the fluctuations transverse to the direction of J grow exponentially till they are 
shut-off by the non-linear back-reaction. The expectation value of the field tends to reach 
the absolute minimum which is on the other side of a potential barrier for values \J\ smaller 
than ^ (Impyf). See Fig. ]!]. 

There are two regimes for the dynamics separated by a critical value of the external field 
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For \J\ > Jc the system quickly reaches a quasiperiodic regime where the field expectation 
value oscillates near the true quantum vacuum. See Figs. |, |, | and ^. 

For \J'\ < Jc the field cannot reach the absolute minimum moving in the longitudinal 
direction except by tunnel effect which is suppressed in the large N limit. Instead, the 
field turns around the potential maximum developing large transverse quantum fluctuations 
(9(m^/A). It takes a time tg ~ ^^^^j^log-^ (spinodal time) for the quantum fluctuations to 
become of the order (9(m^/A). By this time the fleld expectation value makes a spectacular 
jump and starts to oscillate near the true vacuum. See Figs. and |[ 

For later times the evolution in both cases becomes quasiperiodic showing a clear separa- 
tion between fast and slow variables. The fast time dependence of the fleld expectation value 
is explicitly found in terms of elliptic Jacobi functions. Comparison with the full numerical 
solution shows that their parameters exhibit slow time dependence. The smaller is \J'\ the 
slower are the slow variables. We also see that the mass squared is of order \J'\. For \J'\ = 
we recover the known zero external fleld results |^. 

We consider small values of the coupling constant A <^ 1 in order to separate the different 
time scales of the evolution. 



2 Classical field dynamics 

We will flrst propose and solve the problem in classical fleld theory, in order to present results 
that will be used later for comparison. We consider the classical 0(A^)-invariant scalar fleld 
model with quartic self-interaction in the presence of an homogeneous external fleld JT. 
We study in this model the evolution of the ground state after the change J' — J in the 
external fleld. 



2.1 Classical evolution equation and initial conditions 

The action and Lagrangian density are given by. 



S 

c 



d X C 



d^,^{x] 



2 8iV V 



'V. 



We restrict ourselves to a translationally invariant situation. In this case the fleld $ and the 
external fleld J are independent of the spatial coordinates x and only depend on time. The 
evolution equation is then 



and the classical potential 



A 



(2) 



(3) 
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We consider the broken symmetry case (m^ < 0). 

— * 

Choosing the first axis of the internal space in the direction of the external field J ^ 
their components are J = {\/N J, 0, . . . , 0). Consequently, the $ field can be expressed as 

Thus, the potential becomes: 



Kr(0, 07r) 



1 2 r2 A 
- m 0^ H 

2 ^ 8A^ 



(4) 



We choose the ground state as initial state, therefore it has the form (this can be shown 
from the previous expressions) 

(Piti) = 00 ; ^iU) = 



Mu) = ; Mu) = 



(5) 



where 0o is the global minimum of Vj(0). 
The evolution equations for are 



"^' + ^(0'W + 0'^^'W 



0"^.^ = 



(6) 



and with the previous initial conditions, independently of the value of J, the solution is 
simply 

0"^.^ = . (7) 

Therefore, the classical problem reduces to the following one component field problem. The 
evolution equation for reduces to 



0W + 



m2 + ^02(t) 



= +J for t < . 



As it is obvious, the ground state 0o is a constant solution of this equation. We see from 
Eq. @ that 0o, the global minimum of Vj(0o) fulfills Eq. (||). 



We introduce the dimensionless variables 

fx J 



T = \m\ t 



V 



IX 
2 \m\ 



J 



2 \m\ 



Vr 



<0 



A 
2^ 



Vj 



X 1 



2iV ImP 



(9) 
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Therefore, 



since we choose < 0. 

Thus, the classical vacua are the solutions of the cubic equation V^^q^t]) = 0, 



V - ^7 - J = 



(10) 



For small external fields j < 



27 



0.3849002 . . . Eq. (^) has three real roots. The ground 



state (the global minimum) of the potential Vr<o is at 



^ = ^0 = ^ cos 



-arccos 



27j 



1 + i^- _ 3 -2 1 .3 _ 105 .4 5 

2 8 2 128 ^ ^-^ ^ ' ^ ^ 



there is a local minimum (false vacuum) at 



Vfv = ^ cos 



-arccos 



27j 



27r 



and a local maximum at 



2^^ 



105 
128 



= COS 



-arccos 



27j 



+ 



(12) 



For larger external fields j > there is only one real minimum, (the absolute minimum) 



at 



V = Vo = 



1/3 



1/3 



1/3- 



27 



27 f 



3 ^-1/3 81 j5/3 



+ 0{3 



-7/3N 



We shall concentrate in the more interesting case j < 



The field stays in the ground state ?7o with j > from the (negative) initial time Tj till 
r = 0. At r = we suddenly flip the sign of the external field j. 
That is, the time dependence of the external field is 



J 



(j, 0, . . . , 0) for r < 
(-J, 0, . . . , 0) for r > 



(13) 



This change of sign of j introduces some amount of energy in the system. 

As previously stated this problem reduces to a one component field problem. Therefore, 
the potential is given in dimensionless variables by (see Fig. H, P) 



V{ri) 



K>o(r7) 



\ + i^'' ~ iV for r < 
— I rf + \rj'^ + jrj for r > 



(14) 
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The equation of motion for r > is 

^ + (-1 + rf) r] = -J , (15) 

where — 1 + plays the role of an effective mass. The initial state at r = is the ground 
state for r < 0, and it has the form 

v{o) = vo ; m = o. (16) 



2.2 The classical dynamics 

We have shown above that the initial state (the ground state for an external field j) is 
a stationary state of the evolution equations for r < 0. The change j — — j at r = 
breaks this stationarity and the system now finds itself near a metastable state (see Fig. 
H, H). We consider j < when there is a potential barrier (see Fig. |ID. This change in 
the sign of the external field j, increases the energy density of the system by an amount 
K>o('7o) ~ K<o(^o) = 2j?7o- However, even if the field is near the metastable minimum, 
it could have now enough energy to jump above the barrier and reach the global minimum 
of the potential. As the height of the barrier is K->o('7m') the condition to jump above the 
barrier is 

K>o(^?o) > K>o(w) , (17) 

where riu' is the local maximum of the potential K->o(^)- 

Solving this equation one obtains that the field jumps above the barrier provided | j| > jc 
with 



132 + 1575) 

^ = 0.243019.... (18) 

\ IQ'^ 

[Notice that jc < and therefore there is indeed a barrier for these values j < jc-] 
This critical value jc of the external field separates the dynamics in two regimes 

• IjI > Jc- The field jumps above the barrier and reaches the region of the global 
minimum. After this the field continues to jump above the barrier back and forth with 
periodic oscillations. See Fig. |] and dashed line in Fig. |^. 

• IjI < Jc- The system does not jump above the barrier, and oscillates periodically 
between the turning points of the motion 770 and rit around the metastable minimum. 
?7o is given by Eq. ( |Tl| ) and r]t is the solution of the equation 



K>o('7t) = K>o('7o) 

we find, 

. 3 . 11 .2 7 .0 1049 .4 ^, 
r/i = 1 - -J - yj' - -j' - -^3' + O(j^) (19) 

Since both r]t and 770 are positive, r] is always positive. See Fig. ^ and dashed line in 
Fig. 0. 
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3 Quantum field dynamics 



In this section we derive the quantum equation of motion and compute the quantum field 
dynamics using the large N limit. As in the previous section, we consider the 0(A^)-invariant 
scalar field model with quartic self-interaction in the presence of an homogeneous external 
field J'. We study in this model the evolution of the ground state after the change J' — J 
in the external field. 



3.1 Quantum equations of motion and initial conditions 

The action and Lagrangian density are, 



8N V J 



(20) 



where $ is now a quantum operator. 

We restrict ourselves to a translationally invariant situation. In this case the order 
parameter ($(x, t)) and the external field J'{x,t) are independent of the spatial coordinates 
X and only depend on time. Moreover, we choose the direction of J^{t) independent of time. 
Then, we can choose this direction as the first axis in the A^-dimensional internal space. 



J 




0) 



0) 



for t < 
for t > 



(21) 



The field expectation value in the problem we are studying is parallel to the external field: 

$(x) = ia{x), 7f (x)) = i^/N(j){t) + x{x),n{x)) (22) 

with 0(t) = (cr(x)); thus, (x(a;)) = 0. While in the N — 1 directions transversal to the 
expectation value we have (7r(x, t)) = 0. 

The derivation of the equations of motion in the large limit in this case, is analogous 
to the case without external field, that has been explained in detail in [^. Thus, we will 
only present here the main concepts and equations. 

We have one direction parallel to the expectation value, and A^ — 1 transversal directions. 
The transverse fiuctuations dominate in the large A^ limit while the longitudinal fiuctuations 
only contribute to the 1/A^ corrections to the equations of motion. 

The evolution equation for the expectation value is for t > 

{^\x))- 



m H 



N 



-J 



(23) 



where the dot denotes the time derivative. We see in the previous equation that 



X 



{n\x)) 



N 
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plays the role of an effective mass. The last term in Ai'^{t) has a quantum origin and it is 
absent in the classical effective mass. It can be interpreted as an in-medium effect due to 
the presence of vr particles. 

In the Heisenberg picture we can write 

^(^' t)= f -7=7^. Mt) e^'-' + 4 e-''-'] (24) 



V2(27r)= 

where , are the time independent annihilation and creation operators with the usual 
canonical commutation relations. Thus, the (pk{t) are the mode functions of the field and we 
can express (vf ^(x)) as 



(7r2(x)) 1 f d^k 



N 2 J (27r) 



\Mi)\^-Sa 



= l-'-^Mlit) (25) 

where k is an arbitrary renormalization scale, and we will choose k = \mFi\ for simplicity. 
(For details on the renormalization procedure, that leads to the substraction Sd see Ref. [§|.) 

In the previous expressions we have written k, instead of k, since we choose spherically 
symmetric initial conditions. 

The mode functions have the following evolution equations 

ipk{t)+ujl{t)Mt)=0 , ujl{t) = e + Ml{t) . (26) 



We consider values of the field \ J\ < (I'^Pyf) where the potential Vj has two local 
minima. As initial state we take the ground state for t < 

0(0) = 00 ; 0(0) = 

(^,(0) = -^^= ■ <^fc(0) = 7^^. (27) 
/^fc(0) 



where 0o is the global minimum of Vj(0) [defined in Eqs. (P) and (|ri|)] and co'fc(O) 



^2 _|_ _|_ (^_x/2)0o (the square root is real for all /c, because + (A/2)0q > 0). These 
initial conditions correspond to the quantum vacuum (no excitations). 

It is convenient to introduce the following dimensionless quantities [in addition to those 
inEq. (i], 

k A KA 



^,{T) = ^\m\^u{t) ; (7S(r) = — (7r^(t)) ; (28) 
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where m and A stand for the renormahzed mass of the inflaton and the renormahzed self- 
couphng, respectively 0, ||. 

The quantum evolution equation in dimensionless variables are then for r > 0, 



ri{T) + M\t) ri{r) = -J ; (^,(r) + ^J(r) (^,(r) = 



(29) 



with 



M\T) = -l+r]\T)+gJ:{T) ■ 
gi:{r) = g J dq [|^^(r)p - 5(r 
And the initial conditions are 

^(0) = vo ; 

UJq{0) 



q^ + M^r) 



1) M'(t) 



2g3 



(30) 



r]{0) = 
; 0q(O) = -t \j^q{^) 



(31) 



where r^o is given by Eq. ([TT]); and ojqifi) = \J q^ — 1 + rjl because the initial conditions imply 
(7^(0) = 0{g) <C 1. {ujq is real for all g, because rjl > 1). 

It is important to note that the quantum effective mass equals the classical effective mass 
plus an additional quantum term, 5fS(r). This term comes from the in-medium effects due 
to the presence of tt particles. 

In this section we have exposed the quantum equations of motion and the initial condi- 
tions using the Heisenberg picture. The equations shown for the large N limit are the same 
as those obtained in the Schrodinguer picture with a gaussian ansatz for the wave functional 
(as it has been shown in Refs. [151, H) 



^1 



vr; T 



AT^ IT' 



(32) 



where Af^f, is the normalization factor. The relation between the two pictures is made explicit 
by using the relation Ag = —i(f*/ip* . 

gT^ij) gives us the spreading of this wave functional, as we see from Eq. (|28|) (the 
expectation values are independent of the picture). 



3.2 Early time quantum evolution 

At early times, 5'S(r) is negligible in the t] evolution equation because gT,{Q) = 0{g) ^ 1. 
Thus, for early times the quantum dynamics of the expectation value //(r) is the same as the 
dynamics of the classical field presented in section 0. This has been explicitly verified from 
the numerical solution of the full quantum equations 



^-(pl|). As in the classical case, we 
can distinguish two dynamical regimes depending on the value of j compared with jc- 
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Later on, this picture is modified due to the effects of the spinodal instabihties and 
parametric resonances. They make the quantum modes to grow, and therefore so does 
5fS(r). Parametric resonances arrive because A1^(r) oscillates in time. In addition, the 
modes grow exponentially (spinodal instabilities) in the time intervals where 7Vl^(r) < 0. 

These effects have a quantum nature, and they change the dynamics of the expectation 
value ?7(t), due to its coupling with the quantum modes through 5'S(r). 5'S(r) gives the 
spreading of the field wave functional in the transverse directions as shown by Eq. ( p8|) 
(recall that the transverse fields have zero expectation value). 

We comment below on the effects of this phenomena in the dynamics. 



3.2.1 j>j. 

For early times the dynamics of the field expectation value is the same as in the classical 
case. Thus, the system in this case has enough energy to jump above the barrier and reach 
the region of the global minimum of the potential (see Fig. The system continues to jump 
above the barrier back and forth until the effects of the quanta created through spinodal 
instability and parametric resonance become important, i.e., g^ir) ~ 1 (see Figs. 0-§|). 

The system is spinodally instable when A^^(r) < 0, i.e., for //^(r) + 5'S(t) < 1, that for 
early times {gT, ~ 0) reduces to |?7(r)| < 1. It also has parametric resonances due to the 
oscillations of A^^(r). Both mechanisms lead to an abundant creation of quanta transferring 
energy from the 1] degree of freedom to the quantum fluctuations (pk{T)- Therefore, fi'S(r) 
increases while the amphtude of the oscillations of rj^r) decreases. See Figs. ^ and ||. 



3.2.2 j<j. 

The early dynamics is the same as in the classical treatment. For j < jc the system does 
not have enough energy to jump above the barrier and reach the global minimum, and it 
oscillates around the metastable minimum (false vacuum). (See Fig. H). 

The field could go classically to the global minimum region moving in the transverse di- 
rections. But this possibility is excluded due to the initial conditions Eq. (|^): the transverse 
field components and its derivatives are zero and they remain zero by the classical evolution. 
See section 0. 

In the quantum treatment there are in general other two ways of reaching the global min- 
imum region: one is tunneling through the barrier and the second is allowing to circumvent 
the maxima by increasing the quantum probability of finding large values of the transverse 
field components. The tunneling is suppressed in the large limit because it is related 
to the quantum effects in the longitudinal directions which are subleading for large N (see 



comments in section |3.1|) . In the large N limit the leading quantum effects come from the 
transversal field components and those are which we describe. 

In the early (classical like) dynamics, ri{T) oscillates between tjq and rjt given by Eq. ( |T9D 
(see Fig. |^. Thus, as 5'S(t) = the effective mass oscillates between the values (see Figs. 
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A^L. = -l + v'o=J-lf + lf-f + 0{f)>0 (33) 

A^Ln = -l + r/,^ = -3j-if -f/-4/ + O(/)<0 (34) 

Therefore, the low momentum modes grow exponentially in the time intervals where the 
effective mass is negative. This spinodal instability increases 5'S(t) yielding a non-zero 
probability for the transverse field components. 

These instabilities are well approximated by considering a constant negative squared mass 
— /i^ given by the average 

_ ^2 ^ ML^ + -ML. ^ _ 1^.2 _ 9^.3 _ 5^.4 ^ o{f^ ^35) 

This approximation reproduces the growth of 5'S(r) in good agreement with the full numer- 
ical solution of Eqs. (p9D-(pl|) (see appendix |6.1| ) 

^S(r) ^ . (36) 

Thus the probability of finding large values for the transverse field components increases 
with time. 

This solution holds till a time (spinodal time) when gT^^Tg) ~ /x^, this condition gives 
To as the solution of 



~ 77- log i-^)+j-^og{fiTs) = log ( — ^ )+-^log 



'log logi 
logi 



(37) 

After Tg the effect of 5'S(r) in the evolution equations become crucial. 

Soon after t = Tg the order parameter r]{T) makes a spectacular jump and starts to 
oscillate near the global minimum [see Figs. ^, 0. This transition is explained by the fact 
that 5'S(r) = (^^(^))-R ~ 1 implies transverse field components of the order \m\/\/\. 
These large transverse field fluctuations allow the system to turn around the maximum at 
$ = {\/N(j), 07r) = (l"^l\/f Vm' , 0) and reach the region around the global minimum where 
it oscillates quasi-periodically (see Figs. [7|-[Tll). 

The growth of g'^ir) allows ri{T) to reach the region around the global minimum. That 
is, the growth of the transverse quantum fluctuations elliminate the barrier that keeps t]{t) 
near the metastable minimun. 



3.3 Intermediate time quantum evolution 

After the early period described in the previous section, the system enters in a quasi-periodic 
regime. Recall that for different initial conditions the oscillations damped much faster 
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The quasi-periodic regime behaviour found for rj^r) and 5'S(t) suggests that these quan- 
tities are approximately governed by an effective hamiltonian with a few degrees of freedom. 
Actually, we observe from the full numerical solution of Eqs. (p9|)-(pID that 5'S(t) and 7]{t) 
turn to be approximately related as 

gEir) = 1 + a J - (1 - (5 jMt) + jf (38) 

where a and (3 are positive numbers of the order j° and for small g. The coefficients a 
and P are obtained by fitting to the numerical solution. See Table |l|. 



j 


0.05 


0.10 


0.15 


0.20 


0.24 


0.25 


0.30 


a 


0.74 


0.67 


0.60 


0.56 


0.35 


0.29 


0.18 


/3 


1.36 


0.91 


0.60 


0.40 


0.42 


0.49 


0.43 



Table 1: a and /5 values obtained fitting gT,{ri) in the time interval r G [800, 1000] for 
g = 10~^ and various j values. 



We thus find for the effective squared mass, 

M\r]) = -1 + + gj:{r]) =j[a-j+ fif - 2(1 - (3j) 7] + Pv\ 
The evolution equations (|29|)-(|30D in this approximation take then the form, 

i] + M'^{r]) 1] = —j 

We find integrating on 77, 



where, 



Vintiv) 
Eint 



+3 Vintiv) =jE^nt 



Vintim) 



(39) 

(40) 
(41) 

(42) 
(43) 



Notice that Eint depends on the initial conditions, and that rji is a turning point of the 
motion. Eq. (|4l|) can be integrated as follows. 



— rv 
2j{r-n)= / - 



drj 



Eint ^ Vintiv) 



(44) 



with Vi = vi'^i)- 

The fourth order polynomial Eint — Vintiv) has always two real roots rji < V2 corresponding 
to the turning points. Depending on the value of j, the two other roots are a complex 
conjugated pair or two more real roots. For each case: 
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i) a pair of complex conjugate roots rjR ± irji. We then define: 

riR-rii . , _ m 



ii) a pair of real roots rji < rj2 < rj^ < rj^. We then define: 

^2 _ 4:{ri3-r]2){r)4-V3) 



1 



We also introduce two other quantities to simplify the formulae: 



d 



V2 - Vi 



X 



[a - df + 6' 



1/4 



C^Xd2\M\r^,)Tli+j\ 



It is convenient to use in Eq. 



as integration variable u 



The solution of Eq. (4|) can be expressed after calculation as 

{m-vi)[i-cn{C{T-n),k)] 



v{t) =r]i + 



1 + {V2- Vi)X^ - [1 - {V2 - m)X^] cn(C(r - n),k) 
where cn(2;, k) is the Jacobi cosine function, and 



k 



71' 



'1 + 



a — d 



the elliptic modulus. 

The solution (|48| ) oscillates between rji and 772 with period 



T = -K{k) 



(45) 



(46) 



(47) 



(4J 



(49) 



(50) 



where K{k) stands for the complete elliptic integral of the first kind and C is given by Eq. 
(^). The analytical and the numerical solutions are compared in Fig. |ll|. 

The numerical solution would be exactly periodic if Eint were exactly conserved; however 
it is slowly decreasing. The damping is small the smaller is j. See Table |[ 



r 


500 


1000 


1500 


2000 


jEi^tij = 0.05) 


0.02 


0.016 


0.014 


0.011 


J Eint{j= 0.20) 


0.04 


0.025 


0.015 


0.010 



Table 2: E^ for g = 10'^ and j = 0.05 and 0.20. 



Thus, Eint slowly decreases while 77 (r) and Ai'^{r]) oscillate fast showing a clear separation 
between slow and fast variables. We have explicitly found the fast time dependence of the 
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order parameter and the quantum fluctuations in terms of elliptic functions. Comparison 
with the full numerical solution shows that Eqs. (|3^) and (^81) actually reproduce the full 
dynamics provided the constant parameters in Eq. (^) become slow functions of time. 
For example, the parameters displayed in Table D change by 3 — 4% between r = 500 and 
r = 1000. The study of such slow dynamics is a very interesting problem beyond the scope 
of this article and that can be solved using Whitham methods ||10[- 

As is known, for j = the dynamics is governed by a exactly fiat effective potential for 
|?7| < 1 0. Since the slow variables become here slower for j — > 0, they may be governed by 
an effective potential becoming fiat at j = 0. 

Now, we can obtain analytic expressions for 5'S(r) and A4^{t) from the relations g^{f]) 
[Eq. dD] and M^{r]) [Eq. (|3|)], just using the analytical solution r]{T) [Eq. (H)]. The 
relation gT,{ri) indicates that gT, and rj oscillate with opposite phase, up to order j terms; 
gTj{ri) = 1 — rf + 0{j). This 0{j) terms are very important because the terms cancel in 
the effective mass, and the j terms make the mass to be different from zero. This has the 
interesting result of making the time averaged squared mass positive. Such average tends 
asymptotically to a positive value. [Recall that for j = in the broken symmetry case the 
squared mass goes to zero for initial energies below the potential at the local maximum 

(^ = 0) 0,1-] 

The effective mass vanishes asymptotically for j 0. This can be seen explicitly from 
Eq. dD; it implies M'^{t = oo) = 0{j) for small j. 

Notice also that in the j = case the effective mass vanishes asymptotically 0] giving 
the constraint 

1 = + (yfS for j = ; 



and the thick curve in Fig. ^is 0{j) away from such a parabola; showing that the squared 
mass is of order j. 

The mass of the scalars around the ground state m follows from the effective potential 
[see Appendix B]. For j = one finds that the mass rh vanishes ||^, in accordance with 
the presence of — 1 Goldstone bosons. For j ^ the mass becomes non-zero and for 
small j we find in Appendix B, 

= j + 0{j'^ logj) for the ground state. 

The asymptotic mass squared A^^(oo) found above does not coincide with although 
both masses are of the same order of magnitude 0{j). Notice that describes excitations 
around the ground state, while A4'^{oo) corresponds to excitations around states with finite 
energy density above the ground state. 



3.4 Constraints on the trajectory due to the conservation of the 
energy 

As the states considered here are homogeneous and isotropic, the system has an energy- 
momentum tensor with the ideal fluid form. Thus, we can define a dimensionless energy 
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as 



+ J + U^"dq\^,\' + yq'dqq'\^,\'. (51) 

for r > 0. 

The conservation of the energy (for constant external fields) gives us the following con- 
straint in the [rj, gH) plane 

eo > K//d,n(r],S) = -y + J + i^r,'9^ - ^ + ^ + ^ + • (52) 

where eo = V^<o(^o) + 2j^o = K->o(^o) is the energy after the sign change in the external 
field j. 

This constraint is represented as a dotted line in Fig. 0. We see there how the energy 
conservation constrains the trajectories in the (77, (7E) plane. 

Veff dyniv^^) ^an bc interpreted as a dynamical effective potential for the longitudinal 
field 7]. Notice that the equation of motion (^) for //(r) can be written as 

d 

^^^^ ^ -^K//dyn(r?,S) . 

Veff dyniVy ^) defined by eq. (|52D is plotted for j = 0.20 in Fig. |12. 



4 Comments and Conclusions 

Here, we have studied effects of external fields in non-perturbative quantum field dynamics in 
the large limit. The effects of uniform external fields in $^ theory with broken symmetry 
were recently studied in a different framework through classical evolution with random initial 
conditions |jTT| . 



We have studied the evolution of the ground state in an external uniform field J after 
flipping its sign: J ^ — J. We have considered the broken symmetry case and small external 
fields (then the classical potential presents two local minima). 

The change of sign in the external field leads to spinodal instabilities and parametric 
resonances yielding abundant particle creation of the order 1/A. As a consequence, the wave 
functional spreads allowing large values for field components orthogonal to JT. Thanks to 
such large transverse field fluctuations the system goes around the maximum of the potential 
and reaches the global minimum without tunneling. 

The large transverse fluctuations S(r) change the effective mass squared in the evolution 
of the longitudinal field rj^r) [see eqs.(^)-(0)]. In this way rj^r) overcomes the classical 
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barrier and reaches the absolute minimun. The effective dynamical potential Veff dyn{fli ^) 
introduced in eg . (|5^) helps to understand this phenomenon. As shown in fig. 12, the barrier 
in the 77 direction between the two minima of Veff dyniv^ S = 0) dissapears in Veff dyniv^ ^) 
for sufficiently large S. 

After that, the system oscillates reaching an almost periodic regime. That is, fast vari- 
ables (whose time evolution we explicitly solve) oscillate periodically while other variables 
change slowly due to a tiny damping. 

In this quasi-periodic regime the effective squared mass is of order J. Thus, in the J —>■ 
limit we recover a zero effective squared mass consistently with the known results for J = 
where out of equilibrium Goldstone bosons appear [0, ||. 

The damping of the fast oscillations is here much slower than for the oscillations found 
in refs. |^ for other initial conditions. There, the oscillations of A4'^{t) — 7Vl^(oo) were 
damped as 1/t, while the damping is significatively slower for the quasi-periodic solutions 
found here. 
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6 Appendix 

6.1 Appendix A: Spinodal instability 

Here, we obtain the early time solution for the modes in the spinodally resonant band, and 
we estimate the spinodal time, r^. 

Before entering on the calculation of r^, let us recall that the the contribution of the 
spinodal band can be estimated using an average mass for the evolution. The numerical 
calculations show that this estimation is correct, until such contribution becomes of order 
one. 

Hence, an approximate equation for the modes reads. 



The modes with q in the interval between and n are spinodally resonant. There are no 
particles in the spinodally resonant band. Therefore the initial condition for modes in the 
resonant band are: 




(53) 



1 



{q' + \M\0)\) 



1/4 



0,(0) = ~iVQ, = -I (q^ + |-M2(0)|) 



1/4 



(54) 
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The solution of Eq. ( p3D for these modes is: 

1 



T 



2V/i^-g^(g2 + |>i2(o) 1)1/4 



■ (55) 



We thus obtain for the squared modulus neglecting the exponentially decreasing term 



Mr)\' ^ , + (56) 

4/.2 (l - li) ^q2 + \Mm\ 

The contribution of the spinodal band to (7S(r) is given by, 

S,(r)= r (^)|2. (57) 

JO 

Inserting Eq. (^) into Eq. (^) we obtain an estimation for the spinodal growth of the 
quantum fluctuations. To approximately evaluate this integral, we can make further sim- 
plifications in Eq. (|56|). As g//i < 1 and the contribution of the modes with q ^ fi is 
exponentially suppressed, we can expand in to second order in the exponential and to 
zeroth order in the factor outside the exponential. 



lv^g(i")r ~ — e 1^ . (58) 

4/i2 ,/g2 + \M'^i"^' 



In addition, the integrand has its maximum at g = 0(0. lyu) and < /i < 1. Therefore, we 
can approximate fi^ + |7W^(0)| ~ 2/^2 and Jq'^ + \M'^{0)\ ~ fi. Thus, 



MT)\'^j-e'^^e--'^. (59) 



Then integrating over q and using /q^ t>2(iti exp(— t>2r) '^=°° y7r/(16r^, 

g^ir) ^ . (60) 

The spinodal time Tg is, by definition, the time where the instabilities are shut off for all 
< 5' < /i. This happens when the spinodal modes contribution gT^slr) compensates the 
initial (negative) value of M^jj-{t) = — /i^. 

^?S,(r,)^/i2 (61) 

Therefore Tg is given by the following implicit equation. 



+ -^log(/ir,) . (62) 



The spinodal times given by this equation are in good agreement with the numerical results. 



17 



6.2 Appendix B: Ground state results 

The expectation value of the field and the mass of the scalars can be obtained for the 
ground state using the effective potential 0. One finds in the broken symmetry case 
using dimensionless variables, 

2 . -^2 A , 1 
= 1 + m + log — , 

where m is the mass around the ground state of the — 1 scalars. The field expectation 
value is connected with the external source by 

J 

<^ = — • 

Hence, rh? is a function of j defined by the trascendental equation 



A Tfl^ 1 

^"f + ^" + ^2^^^^2=3- (63) 

For j = one finds that m = in accordance with the presence of iV — 1 Goldstone bosons. 
For j 7^ the mass rh? becomes non-zero and for small j we find 

= J + 0{j'^ logj) for the ground state. 

[Notice that Eg. (|6^) becomes the classical equation Eq.(|IO|) for A = setting rri^ = —l + r]"^, 
as expected.] 
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CD C\J O OvJ CD 

O O O CD CD CD 

A 

Figure 1: Classical potential at r < for small field j = 0.20 < 2/v^27- The potential 
presents two minima (the ground state r/o and a metastable false vacuum state rjfy) separated 
by a potential barrier. 



20 




A 

Figure 2: j = 0.25 > jc- Classical potential at r = and for r > 0. The change of the sign 
of the external field at r = increases in 2jrio the energy density of the system, and the 
system has for j > jc enough energy to jump above the barrier. 
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Figure 6: j = 0.20 < jc- Classical potential at r = and for r > 0. The change of the 
sign of the external field at r = increases in 2j?7o the energy density of the system but the 
system does not have enough energy to jump above the barrier. 
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Figure 12: j = 0.20 < jc- Plot of the dynamical effective potential Vejf dyniv^^) given 
by eg. (|52D . Veffdyniv^^) is bounded from above by the total energy. This constrains the 
trajectory represented in Fig. 0. 
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